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Abstract. This paper discusses recent results and near-term prospects of the long- 
baseline neutrino experiments MINOS, MINOS-I-, T2K and NOvA. The non-zero 
value of the third neutrino mixing angle ^13 allows experimental analysis in a manner 
which explicitly exhibits appearance and disappearance dependencies on additional 
parameters associated with mass-hierarchy, CP violation, and any non-maximal 023- 
These current and near-future experiments begin the era of precision accelerator long- 
baseline measurements and lay the framework within which future experimental results 
will be interpreted. 
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1. Introduction to Long-Baseline Accelerator Experiments 

1.1. Motivation and 3-Flavor Model 

Beginning with the snccessfnl operation of the K2K experiment |23], the physics com- 
mnnity has seen a profonnd expansion of onr knowledge of the mixing of nentrinos, 
driven by long-baseline accelerator experiments [ZmiEllllSlES], experiments stndying 
atmospheric nentrinos pUl (23, [ 2 H] , solar nentrinos [ 2 ] , and, most recently, high-precision 
experiments with reactor neutrinos [55]. In this article we describe the current gener¬ 
ation of the running long-baseline neutrino experiments T2K and MINOS/MINOS-I- , 
and the status of the NOvA experiment which was commissioned in 2014. Each of these 
experiments was designed with primary and secondary goals. For example, MINOS had 
as its principle justihcation the measurement, via disappearance of of the mixing 
parameters sin^ 6'23 and Am| 2 , with particular emphasis on Am‘^ 2 - T2K and NOvA 
were primarily designed to elucidate the structure of the neutrino sector by studies of 
Ue appearance. However they are making and will continue to make very signihcant 
contributions to the study of disappearance as well. Similarly, MINOS has measured 
Ue appearance and the angle 6 * 13 . 

This situation leads us to a point we will emphasize throughout this article, namely 
that the traditional distinction between various modes of study of mixing, matter effects, 
and CP violation is rapidly giving way to a more integrated approach which utilizes both 
major types of signals to gain maximal information about the somewhat complicated 
3-neutrino sector. We will hrst discuss the relevant formalism and the major measure¬ 
ments required to test it. Next we will provide a technical overview of the powerful 
neutrino beams required for the measurements, and then proceed to discuss the experi¬ 
ments together with their current measurements and expected sensitivities. Finally, we 
conclude the article with a discussion of the near-term future for these experiments. 

It is our goal to familiarize the reader with the surprisingly rich information already 
available in studies of this sector of physics, the only one currently not well-handled by 
the Standard Model. In addition, we will provide context for future discussions of 
progress to be provided by these experiments and by the exciting future world of very 
large experiments, and very long baselines. 

1.2. Core measurements 

In the three-neutrino model there is a close relationship among the disappearance and 
appearance modes of oscillation study, going back to their origin in the PMNS matrix. 
Following reference [36], it is possible to write the disappearance possibility for muon 
neutrinos in vacuum, as 
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Am^ L 

= 1 - 4sin^ 023 cos^ 013(1 - sin^ 023 cos^ 0i3) sin^ — (1) 

where Arn^^j incorporates effective leading dependences on the additional PMNS pa¬ 
rameters Am|i, 013 , and 5cp as 


Amlfj: = Am 32 -|- Am^i sin^ 0 i 2 -|- Am^i cos 5cp sin 0 i 3 tan 023 sin 20 i 2 

Equation may be simply manipulated to yield a form appropriate for the baseline 
of the T2K experiment, namely 


/\ ?T7 ^ T 

~ 1 - (cos"^ 013 sin^ 2023 + sin^ 20i3 sin^ 623 ) sin^ (2) 

Here we see the vital role of the mixing angle 0 i 3 , which couples in the PMNS matrix 
to the CP-violating phase 6 cp. It is now well-known that this angle is relatively large, 
approximately 9 degrees [9]. T2K, MIN0S/MIN0S-1-, and future NOvA measurements 
of this angle and its consequences are discussed in sections lllE]and[3 

Because the earth between the beam creation point and the detector location forms 
an essential part of any long-baseline experiment, its effects on the measurements must 
be considered. This creates both problems and opportunities - problems because of the 
introduction of degeneracies between matter effects and CP violation, and opportunities 
because of the possibilities to exploit the differences between neutrino and antineutrino 
interactions, and from the two mass hierarchies. 

1.2.1. Appearance measurements and sire^ 29 13 

A large value of 0 i 3 is key to allowing an integrated approach to oscillation studies. 
MINOS, T2K, and soon NOvA use the appearance channel for z/g with a muon neutrino 
beam to probe 0 i 3 directly. In appropriate approximation for a muon neutrino with 
energy E^, of 0(1) GeV traveling a distance of 0(100) km, the leading order equation 
governing the appearance probability is: 

sin" sin" 2»,3 sin" NAt (3) 

Equation is applicable for both neutrino and anti-neutrino oscillations. 

The difference between neutrinos and anti-neutrinos in oscillation appears as a sub¬ 
leading effect, including the solar parameters 0 i 2 and Am^i and CP violation phase 5cp. 
The probability is expressed |32l E3] as 

P{^fi Ug) ~ sin^ 2013 T 1 — a sin 2613T2 + a sin 2013 T 3 -|- 


(4) 
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where a = is the small (~ 1/30) ratio between the solar and atmospheric squared- 
mass splittings, and 


Ti 

T 2 

T3 

T4 


sin2[(l-A)A] 


sin 623 


(1 - A2) 


sin hcp sin 26^2 sin 2623 sin A 


cos hcp sin 2012 sin 2023 cos A 


. 2 . sin (AA) 


sin(AA) sin[(l — A)A] 
A 1-A ’ 

sin(AA) sin[(l — A)A] 


A 


1 - A 


cos 023 sin 2012 - 


A2 


Here, A = and A = 2^/2GpngEy /where N^, is the electron density of 

Earth’s crust. In Equation the sign of the second term changes for anti-neutrinos, 
governing CP violation when all three mixing angles, including ^ 13 , have nonzero values. 
With current best knowledge of oscillation parameters, the CP violation (sub-leading 
term) can be as large as ~30% of the leading term. 

The A dependence arises from matter effects (caused by additional terms in the 
Hamiltonian for the electron component of the neutrino eigenstate), which are coupled 
with the sign of Am 32 . In this paper, we refer to Am ‘^2 > 0 as the normal mass hierarchy 
and Am ‘^2 < 0 as the inverted one. 


1.2.2. Hierarchy, Octants, and CP Violation (5cp) 

The subleading terms shown in Equations [T] and make it possible in principle 
for oscillation measurements to be sensitive to the octant of 6*23 (023 < 7 ^/ 4 : or 
023 > 7 r/ 4 /, in the case it is not maximal (i.e. = vr/4). Even for maximal mixing, 
the additional dependencies in Equation and further terms added to Equation 
enable the combination of disappearance and appearance results to begin to give clues 
to the CP and hierarchy puzzles. 


2. Beamlines 
2.1. NuMI 

The conceptual beginning of the NuMI beamline dates to the era of the construction 
of the Fermilab Main Injector. The beamline, together with its associated tunnels, 
experimental halls, surface buildings and infrastructure, was built between 1999 and 
2004. Datataking with the beam began in March of 2005, and continues to this time. 
The complex currently consists of a primary beam transport, the NuMI target hall, a 
675 meter He-£lled decay pipe of one meter radius, a hadron absorber and muon flux 
monitor area, meters of rock shielding, and two experimental halls - the first housing the 
MINOS near detector and the MINERvA detector, and the second housing the NOvA 
near detector. Figure shows the general conhguration of the beamline. 
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Figure 1. Layout of the NuMI beamline at Fermilab, showing main components of 
the target, focusing, and decay systems [5]. 


Beam from the Fermilab Booster is accelerated in the Main Injector to 120 GeV, 
and then extracted with a system of fast kickers. As part of an npgrade to beam power, 
commissioning is nnderway to stack beam in the Fermilab Recycler Ring before transfer 
to the Main Injector. At the NnMI target the proton beam consists of 6 batches, with 
a total extraction period of 10 pis. The time between extractions has varied from 2.2 s 
to the cnrrent 1.3 s. This low dnty factor allows the MINOS and NOvA experiments to 
trigger on a simple timing window, which facilitates the snrface location of the NOvA 
far detector. Overlaying mnltiple injections from the booster in the main injector (slip¬ 
stacking) has allowed the beam intensity to reach 375 kW. The NnMI beam for NOvA 
is anticipated to reach 700 kW, with a similar time strnctnre. 

Secondary hadrons created in the interaction of the extracted proton beam with a 
94 cm graphite target are focused by a system of two magnetic horns. Historically, the 
focused hadron momentum (which translates into neutrino energy) has been adjustable 
by moving the target w.r.t the hrst horn. This is more difficult in the NOvA beam 
configuration, which is optimized for the off-axis application. Figure shows the 
focusing schema of the NuMI beam.The majority of the data samples used in MINOS 
physics analyses have used the NuMI beam line focused at its lowest practical energy 
configuration, with a peak neutrino energy at approximately 3 GeV. Additional samples 
at higher neutrino beam energy settings have been used to provide information on the 
intrinsic z/g component of the beam, an important and irreducible background for the 
measurement of sin^ 26 ^ 13 . Overall beam production for the MINOS running period 
2005-2012, is shown in Figure]^ 

An important feature of horn-focused beamlines like the NuMI beamline is the 
ability to convert from focusing positive hadrons (primarily tt’*') to negative hadrons 
(primarily 7 r“), which creates a beam heavily enriched in antineutrinos. Using this 
beam, MINOS has published special studies of the oscillation parameters of |4D] . 

Targetry which can withstand the repeated high power proton pulses needed for 
a neutrino beam represents a technical challenge. The NuMI targets used for MINOS 
data taking were constructed of 47 segmented graphite fins. A total of 7 targets were 
used in the period 2005 to 2012, with exchanges usually due to failures in auxiliary 
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Figure 2. Detail of the magnetic focusing horn system for the NuMI beam line mi- 
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Figure 3. Total NuMI beam delivery during MINOS running. The total collected 
in the principle neutrino configuration was 10.7 x 10^° protons. For the antineutrino 
configuration the total was 3.4 x 10^° protons. 


cooling systems. In one case the target material experienced significant degradation, 
visible in the produced neutrino rates. Significant engineering changes have occurred 
for the targets to be used in the NOvA beam, which must withstand 700 kW operations. 
These include detailed changes to the graphite fins to allow for an increase in primary 
beam spot size from 1.1 mm to 1.3 mm, and, importantly, a significant relocation of the 
water cooling tubing to decrease its vulnerability. With these changes, it is expected 
that the NuMI targets in the NOvA era will survive a minimum of a year of high power 
operation before any replacement is needed. 

2.2. Off axis neutrino beam 

An off axis neutrino beam (OAB) configuration [I] is a method to produce a narrow 
band energy neutrino beam. In the OAB configuration, the axis of the beam optics is 















































































Neutrino Oscillations with the MINOS, MINOS+, T2K, and NOvA Experiments. 7 


intentionally shifted by a few degrees from the detector direction. With a hnite decay 
angle, the neutrino energy becomes almost independent of the parent pion energy due 
to characteristics of the two body decay kinematics of the pion with Lorentz boost. The 
off-axis beam principle can be illustrated with a simple algebraic example. Let us model 
the beam as consisting of pions which are fully focused in the on-axis directions. The 
transverse and perpendicular components of the decay neutrino momentum obey the 
relations: 


Pt = P*sin9*,PL = iP*{I + cos 9*) ^ (5) 

where P* and 9* are the decay momentum and angle in the rest frame of the 
decaying particle. The fixed off-axis angle condition is 9iab = Pt/Pl- 

Near 9* = 7r/2, we have APt ~ 0 for variations in 9*, and therefore 
APt 

APl = ^^0 ( 6 ) 

aiab 

Physically, the constraint on the angle means that the variation of neutrino energy 
that normally occurs when 9* varies is greatly reduced, and parent particles of many 
energies contribute to a single peak in neutrino energy. As a consequence, the peak 
energy of the neutrino beam depends on the off-axis angle. Figure illustrates this 
effect graphically for several off-axis angles in the NuMI configuration. 

By changing the off axis angle, it is possible to tune the neutrino beam energy to 
maximize the sensitivity of the oscillation parameters. As a reference, the off axis angle 
can be varied from 2.5 to 3.0 degree in the T2K beamline, which corresponds to a mean 
energy of neutrinos in the range from 0.5 to 0.9 GeV. 

The neutrino energy spectra at the far detector (Super-Kamiokande) with different 
off-axis angles in T2K are shown in Figure 00- In T2K, the off-axis angle is set to 
2.5 degrees. 

The NOvA experiment is situated at an off-axis angle of 14 mrad (0.8 degrees). 
With a higher beam energy focusing than used for the MINOS program, this results in 
a large flux at the neutrino energy associated with oscillations. At the same time, it 
reduces backgrounds from neutral current (NC) interactions from higher energies and 
from intrinsic beam z/g, which have a wider energy distribution. 


2.3. T2K neutrino beam 

J-PARC, the Japan Proton Accelerator Research Complex, is the accelerator complex 
supplying 30 GeV protons to the T2K experiment. An intense neutrino beam with a 
narrow-band energy spectrum is produced using the off-axis technique. The beam energy 
is tuned to the oscillation maximum (~ 600 MeV for the T2K baseline of 295 km), which 
also suppresses the high energy component contributing to background generation. The 
left plot in Figure shows the prediction of the T2K neutrino beam flux at the far 
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Figure 4. The neutrino oscillation probability of —>■ and the neutrino energy 

spectrum with different off-axis angles in T2K from [3] 



E, (GeV) 


Figure 5. The neutrino spectra (flux times cross-section) for various angles in a 
medium-energy NuMI beam [39]. 
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detector, Super-Kamiokande (Super-K). The flux is dominated by muon neutrinos with 
a small fraction (at the level of a few %) of intrinsic electron neutrinos, referred to as 
’’beam z/g”. The beam z/g component is a major background when searching for electron 
neutrino appearance. 


run 1-4 at SK 



run5c at SK 



Figure 6. Prediction of the T2K beam flux for neutrinos and antineutrinos at Super- 
K. The left plot is for the neutrino beam mode made by focusing the positively charged 
particles, and the right is for the anti-neutrino mode made by focusing the negative 
ones. The flux above E,, = 10 GeV is not shown although the flux is simulated up to 
E,, = 30 GeV. 


The proton beam is directed onto a graphite target which is designed to accept 750 
kW beam power. The target is a graphite rod of 91.4 cm long and 2.6 diameter with 
l.dig/cniA density. The target is helium-cooled. Since the current beam power is still 
around 350 kW, there is still a margin for safety. The details of the T2K target are found 
in ref. [TB]. The positively charged particles (mainly pions) produced are focused by 
three magnetic horns, typically operated at 250 kA. The decay of the charged particles 
in a 100 m decay volume produces the neutrino beam. By reversing the direction of 
the horn current, negatively charged particles are focused to produce the anti-neutrino 
beam. The prediction of the anti-neutrino beam flux is shown at the right plot in 
Figure Thanks to the off-axis technique, the signal to noise ratio of the anti-neutrino 
beam flux is as good as 30 [|] at the flux peak, while the wrong sign component of 
neutrinos is broadly distributed in energy. 

T2K started physics data taking in January 2010. Although the data taking was 
interrupted on March 11, 2011 by the Great East Japan Earthquake, the experiment 
collected 6.57 x 10^° POT (protons on target) for analysis before May 2013 [§| The 

f The signal to noise ratio of anti-neutrinos to neutrinos is typically much worse than that of neutrinos 
to anti-neutrinos. 

§ J-PARG stopped operation in May 2013 because of the hadron hall accident. 
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history of data taking is shown in Figure In 2014, anti-neutrino beam running began. 
Today, a maximum beam power of 370 kW has been recorded in J-PARC. 

- Total Accumulated POT for Physics 

• v-Mode Beam Power 



Figure 7. The history of the delivered protons to the T2K experiment for analysis. 
The dots show the number of protons per pulse, and the lines show the integrated 
number of protons. The red dots are for the neutrino beam running, and the purple 
dots are for the anti-neutrinos. 


3. Electron Appearance Analysis 

3.1. MINOS 

The MINOS experiment, with its magnetized steel calorimeters, was principally designed 
to detect and classify the charged current reactions -b A —>■ -f A and -b A —)■ 
-b X. In order to detect and measure the appearance of z/g and z/g, which indicate 
a non-zero ^ 13 , sophisticated statistical techniques must be used to disentangle the 
relative contributions of this signal from the similar neutral-current background. To do 
this, MINOS uses the LEM (Library Event Matching) technique [H]. This procedure 
uses large (> 10 ^) simulated samples of signal and background events to form event- 
by-event comparisons of the observed deposited charge in detector channels with the 
equivalent simulated deposited charge in the library events. The LEM procedure gives 
a set of output variables which are used as input to a simple artificial neural net, giving 
a statistical discriminant ( a^EM ), which can be used to identify signal and background 
components of the data. 

The discriminant a^EM is formed by the output of a neural net which has been 
given as inputs the reconstructed event energy and characteristics of the 50 best-matched 
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library events, namely i) the fraction of these library events that are nue CC events, 
(ii) their average inelasticity (y) and (iii) the average overlapping fraction of charge on 
strips between the data and the 50 library events. Events with ulem > 0.6 are selected 
for farther analysis. 

Next, to search for appearance of z/g dne to the oscillation phenomenon, the spectra 
of the varying backgronnd components present in the beam (z/^-CC, NC, and residnal 
beam i>f,-CC) need to be estimated nsing data from the near detector. This is done in 
MINOS by comparing samples obtained from different beam focnsing confignrations for 
decaying secondaries, as discussed in ref. |13]. Figure from ref. |15] shows the far 
detector MINOS data and the expected backgrounds, for various bins of olem- 

The final elements required to produce an appearance measurement are extrapolation 
of the background (and oscillated signal) estimates between near and far detectors, and 
an estimate of the signal efficiency. The first is done by comparison of the background 
measured in the ND with its simulated value, giving a correction factor that can be 
applied to the equivalent simulation of the far detector. The technique is simpler than 
that used for MINOS CC appearance measurement; however the essential equivalence 
of the methods has been demonstrated in ref. HU. In order to estimate the signal 
efficiency, hybrid events were created by substitution of a simulated electron shower in 
shower-subtracted, well-identified CC events. The efficiencies obtained were > 55% in 
both beam configurations. 

The principal systematic errors affect the result are uncertainties in the background 
estimation and in the signal efficiency. They are 3.8% (4.8%) and 2.8% (3.1%), re¬ 
spectively, for the u (p) modes. The measurement is dominated by statistical errors, 
affecting both the signal and the background estimation. MINOS systematic errors are 
discussed further in section 15.31 


After establishment of these techniques and their systematic errors, MINOS can now 
use the near detector spectrum to extrapolate the expectation of signal and background 
for hypothesized values of the the physical parameters 6 * 13 , 6 cp, and mass hierarchy, to 
determine statistically allowed and disallowed regions. The overall background estima¬ 
tion for the z/ beam configuration is 127.7 background events. For parameter values of 
sin^ 26 'i 3 = 0 . 1 , 6 cp = 0 , ^23 = vr/4 and normal hierarchy, 33.7 ± 1.9 appearance events 
are expected, giving a signal/background ratio of S/B = 0.26. A total of 152 events 
are observed. Contributions to the analysis from the P beam configuration are small, 
totaling only 21.4 expected and 20 observed events. The result from ref. |1^ is shown 
in Figure for the full MINOS dataset, consisting of 10.6 x 10^° POT for the z/ beam 
conhguration and 3.3 x 10^° POT for the P beam configuration. This analysis does not 
distinguish between neutrinos and antineutrinos. The result shows the characteristic 
features of periodic variation with dcp, as well as shape inversion with hierarchy choice. 
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Figure 8. MINOS far detector data with statistical errors [JS] used for appearance 
analysis, compared with expectations for sin^(26*i3) = 0.051, Ato |2 > 0, 5 = 0, and 
6*23 = 7^/4 


MINOS cites best-fit values for 2 sin^ 2^13 sin^ 6^23 of 0.05llog3Q under the normal hierar¬ 
chy assumption, and 0.093 ]'^qq| 9 for the inverted hierarchy, with 90% conhdence ranges 
of 0.01 - 0.12 and 0.03 - 0.18 , respectively. The best hts are all computed for 5cp = 0, 
and 6*23 < 7 r/ 4 . 


It is of interest to examine the parameter space probed by this appearance analysis 
more closely. The MINOS collaboration has computed the change in likelihood for ex¬ 
cursions of the CP-violating phase dcp for four combinations of the hierarchy and 6*23 
octant parameters. This result from ref. |1S] is shown in Figure 10 
disfavors 31% of the total 3-parameter space {5cp, 


The experiment 


hierarchy, octant) at 68 %C.L. and 
shows a suggestive, but statistically limited preference for the inverted mass hierarchy 


scenario. 
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Figure 9. MINOS allowed contours at 68 % and 90% confidence for the measured 
quantity 2 sin^ 20 i 3 sin^ 6 * 23 , as a function of 5cp. Systematic and statistical 
uncertainties are included, and results for both assumed hierarchies are displayed. 


3.2. T2K 

The first evidence of non-zero ^13 was reported [a m] in the —)■ appearance 

channel in 2011. Today, with more data collected in T2K, the —)■ z/g transition is well 
established [7]. Twenty-eight electron candidate events in T2K have been observed in the 
T2K far detector (Super-K) by requiring one Cherenkov ’’ring”, identified as an electron 
type with visible energy greater than 100 MeV. In addition, a newly developed algorithm 
was applied to suppress background events with a vr^ —)■ 27 , where one of the photons is 
missed in reconstruction. The details of event selection are found in The number 

of observed events, compared with the expectations is shown in Table The observed 
number of events, 28, is significantly larger than the expected number, 4.92 ±0.55, with 
013 = 0 , but is consistent with the expectation of 21.6 with sin^ 20 i 3 = 0.10 and ^cp = 0 . 

The best £t value of 0 i 3 has been evaluated to be sin^ 20 i 3 = 0.140 ± 0.038(0.170 ± 
0.045) with a 68 % conhdence level (C.L.), by fixing the other oscillation parameters: 











Neutrino Oscillations with the MINOS, MINOS+, T2K, and NOvA Experiments. 14 



Figure 10. Variation of likelihood (compared to best fits) with 5cp, plotted as 
—2Aln{L), for the observed Pg appearance in MINOS, shown for varying combinations 
of other oscillation parameters. Values above the horizontal dashed lines are disfavored 
at either 68% or 90% C.L. 


Table 1. The number of observed events with the MC expectations and efficiencies. 
The oscillation parameters are assumed to be sin^20i3 = 0.1, sin^ 623 = 0.5, |Ato| 2 I = 
2.4 X 10“^ eV^, dcp = 0, and Ato §2 > 0 


Data 

Total 

MC 

Signal -)■ z/g 
CC 

CC 

Beam z/g + z/g 
CC 

NC 

z/g events 

28 

21.6 

17.3 

0.1 

3.2 

1.0 

Efficiency [%] 

- 

- 

61.2 

0.0 

19.1 

0.4 


sin"^26'i2 = 0.306, Am 2 i = 7.6 x 10“^ eV^, sin^ 6^23 = 0.5, |Am 32 | = 2.4 x 10“'^ eV^, and 
dcp = 0. Figure 11 shows the electron momentum versus angle distribution (sensitive 
to the oscillation), which is used to extract the oscillation parameters sin^ 26^13 and (5cp 
which give the best £t values. The signihcance for a nonzero 6*13 is calculated to be 


7.3 a. 

Allowed regions for sin^ 2^13 as a function of dcp are evaluated as shown in Figure [T^ 
where the values of sin^ 623 and Am ‘^2 are varied in the £t with additional constraints 

















Neutrino Oscillations with the MINOS, MINOS+, T2K, and NOvA Experiments. 15 



5 0 500 1000 1500 

Momentum (MeV/c) 


Figure 11. The T2K electron momentum versus angle distribution for 28 single¬ 
ring electron events, together with the MC expectation in [7]. The best fit value of 
sin^ 2013 = 0.140 in the normal hierarchy case is used for the expectation. 


from [H]. In order to be sensitive to (5cp, T2K uses the value of 6 * 13 , 0.098 ± 0.013, from 
reactor experiments in PDG2012 [9]. The —2AlnL in the £t as a function of 5cp is 
extracted, and is shown in FigureThe T2K measurement, together with the reactor 
013 value prefers hcp = — 7 r /2 with an exclusion of 0.197r < hcp < 0.807r(—vr < hcp < 
— 0 . 977 r and — O.Odvr < hcp < tt) with normal (inverted) hierarchy at 90 % C.L. This 
may be a hint of CP violation in neutrinos. 

4. Precise measurements of oscillation parameters 

4 . 1 . T2K 

The most precise measurement of 023 has been carried out [1] by T2K based on the data 
set of 6.57 X 10^° POT. First, a single-ring muon sampl^ is selected by requiring one 
muon-type Cherenkov ring with momentum greater than 200 MeV/c in Super-K. The 
details of the event selection are found in |3]. One hundred twenty events are selected 
while the expectation without neutrino oscillation is 446.0 ± 22.5 (syst.). The neutrino 
energy for each event is calculated under the quasi-elastic (QE) assumption using the 

II The number of rings corresponds to the number of observed particles in Super-K 
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sin^20j3 


Figure 12. The allowed regions for sin^ 20i3 as a function of (5cp for the normal mass 
hierarchy case (top) and the inverted one (bottom) from [7]. The value of 0i3 from 
reactor experiments in PDG2012 is shown as the shaded region 


expression 

^rec ^ ~ ~ + 2(m„ - Eb)Ei 

2{mn - Eb - El + Pi cos 6 i) 

where rup is the proton mass, the neutron mass, mi the lepton mass, Ei the lepton 
energy, and Eb = 27 MeV the binding energy of a nucleon inside a nucleus. Figure [l^ 
shows the neutrino energy of the observed 120 events with the MC expectations for 
neutrino oscillations. 

Using the number of events and the neutrino energy spectrum, the oscillation 
parameters (sin^ 6 ^ 23 , Am 32 ( 33 )) are estimated with an un-binned maximum likelihood 
£t for the normal (inverted) mass hierarchies. Details of the method are found in |1]. 
The result is shown in Figure [Thj 

The best fit value with the ID 68 % conhdence intervals are sin^ 6*23 = 0.514((Q;g5g 
(0.511 ± 0.055) and = 2-51 ± 0.10 = 2.48 ± 0.10) x lO'^eV^ for the 

normal (inverted) hierarchy case. The result is consistent with the maximal possible 
disappearance probability and is more precise than previous measurements, especially 
for sin^ 623 . 
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Figure 13. The — 2A In L value as a function of 5 cp with the reactor 6*13 constraints 
for normal and inverted hierarchies from [7]. The likelihood is marginalized over ^ 13 , 
023 and Am| 2 - The 90 % C.L. is evaluated by using the Feldman-Cousins method. 
The (5cp regions with values above the 90 % C.L. lines are excluded. 



Reconstructed v Energy (GeV) 


Figure 14. The neutrino energy spectrum for single-ring muon events in [7]. (Top) 
The observed and expected spectra with the event categories in the simulation. 
(Bottom) The ratio of the observed spectrum to the no-oscillation hypothesis with 
the best-fit case. 
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Figure 15. T2K Contours of oscillation parameters sin^ 023 versus | Ato 32 (^ 3 ) | for 68 % 
and 90 % C.L. region. The one-dimensional profile likelihoods are also shown for each 
oscillation parameter in the top and right windows. The ID —2A In ^critical values for 
the normal mass hierarchy are shown in the windows. The Super-Kamiokande m and 
MINOS [S3] results are also shown for comparison. 


4-2. Joint analysis of Un disappearance and z/g appearance samples in MINOS 


Beginning in 2005 (2003 for collection of atmospheric data), the MINOS experiment 
has provided precision measurements of the oscillation parameters Am^ and sin^(26') for 
effective dehnitions of these parameters in a two-neutrino approximation. Most recently, 
MINOS has quoted 2.28 x 10“^ < |Am 32 | < 2.46 x 10“^ eV^ ( 68 % conhdence) and a 
90% C.L. range for 023 of 0.37 < sin^ ^13 < 0.64 (both normal mass hierarchy), using 
a complete 3-neutrino description of the data [5lj. Atmospheric data and appearance 
data are included in the combined hts. In particular, the z/g appearance data and 
atmospheric data provide, in principle, sensitivity to additional information concerning 
mass hierarchy and CP phase. As an example, the atmospheric data sample, divided into 
neutrino and antineutrino samples for up-going multi-GeV events, shows different matter 
effects for normal and inverted hierarchies. In the current sample, these additional 


sensitivities are limited, as shown by the presentation in Figure 1^, and the htted 
values of | Am‘^ 2 1 sin^ 623 are consistent 


The £t results are obtained using 
constraints from external data. In particular, a value of sin^ ^13 = 0.0242 ± 0.0025 
has been taken from a weighted average of reactor experiment values |55], and solar 
oscillation parameters are taken from 
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Figure 16. MINOS 3-flavor oscillation parameters using both disappearance and 
appearance data. Left: Confidence level contours for assumed normal and inverted 
hierarchies, computed using —2Aln{L) w.r.t. the overall best fit point (star). Right: 
One-dimensional likelihood profiles for the parameters. All results from |54) . 


4-3. Joint analysis of Un disappearance and v^. appearance samples in T2K 


The oscillation probability of —?• Vf, depends on many oscillation parameters: sin^ ^ 13 , 

sin^ 6 * 23 , ^^32 ^cp; fhat of —?• depends mainly on sin^ 6*23 and Aml 2 - Therefore, 

all oscillation parameters can be efficiently extracted by htting two data samples of 
and z/^ simultaneously. For this purpose, the T2K collaboration developed analysis 
techniques to fit both z/g and z/^ samples. One method is based on the A log L method, 
and the other is based on the Markov Chain Monte Carlo (MCMC) method. The 6*13 
constraint from PDG2012 [9] is applied in the analysis. 

With the AlogL method, T2K measures sin^ 26 * 13 , sin^ 6 * 23 , and 5 cp as 

shown in Figure [T7| The results are consistent with those shown in Section 3.2 and 4.1 
and the correlations between parameters are properly treated. 

With the MCMC method, the quantities —2A In L{= In L((5cp)~ln T(best fit values)) 
in the £t, as a function of 5cp, are evaluated as shown in Figure f^ The best £t value 
and the preferred regions at 90 % C.L. are consistent with the result of the Ve only 
sample shown in Figure [T^ 

In Figure [T^ the credible intervals calculated in the MCMC method are shown in 
the sin^ 6*23 versus |Am 32 | plane, for both normal and inverted mass hierarchy cases. 
The results are compared with other measurements by Super-Kamiokande na and 
MINOS [5l]. The T2K best fit point is found to lie in the normal mass hierarchy as 
shown in Figure [T^ 
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Figure 17. T2K Contours of oscillation parameters calculated with the AlogL 
method for both normal and inverted hierarchy cases [33]: (top-left) sin^ 023 versus 
|Ato 32 (i 3 )|, (top-right) sin^ 0i3 versus |Am 32 ( 33 )|, (bottom-left) sin^ 0 i 3 versus (5cp, 
and (bottom-right) sin^ 023 versus sin^0i3. The 90 % (68 %) C.L. are shown in the 
solid (dashed) lines with the best ht point shown by the mark. 


5. Systematic Uncertainties 
5.1. T2K Beam 

In accelerator neutrino beam experiments, understanding of the properties of the 
neutrino beam is very important. An experiment is usually designed to cancel hrst- 
order uncertainties of the neutrino beam by adopting the ’’two detectors” technique, in 
which one detector, located near the beam production point, is used to monitor the beam 
and the other, far detector, studies neutrino oscillations. By normalizing the neutrino 
events with the near detector measurement, the systematic uncertainties of the beam 
and the neutrino cross sections are largely canceled. Even with the cancellation, for a 
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Figure 18. (Left) The T2K credible intervals in the sin^ 6*13 versus 5 cp plane 
calculated in the MCMC method [34] . The 68 % intervals are show by the dashed lines 
and the 90 % intervals are by the solid lines. These are constructed by marginalizing 
over both mass hierarchies. The red lines are extracted only from the T2K data set, 
and the black line and the point are extracted from both the T2K data set and the 
reactor @13 constraint. (Right) The —2A InL value as a function of 5 cp in the MCMC 
method with both and samples and reactor ^13 constraint [34] . The Jcp regions 
with values above the 90 % C.L. lines are excluded. 


precision measurement, understanding of the beam itself is essential. 

In the T2K experiment, the neutrino beam is simulated by incorporating real 
measurements of hadron production, among which large contributions come from the 
CERN NA61 experiment [101 E]. The uncertainties of hadron production decaying 
into neutrinos are directly related to the systematic uncertainty in the neutrino beam. 
The uncertainty of the neutrino beam flux at the far detector in the T2K experiment 
is shown in Figure 20 The fractional uncertainty of the beam is at the 10 ~ 15 % 


level, of which the largest component is still the uncertainty in hadronic interactions. 
The second largest component is due to the combined uncertainties of proton beam 
parameters, alignment of the beam line components and off-axis angle. In T2K, the 
beam stability and the off-axis angle are directly monitored using the neutrino beam 


monitor: INGRID [T2]. As shown in Section 5.2, the relative beam flux uncertainty 
between the near and far detectors can be reduced to the 3 % level including the 
uncertainty in cross sections constrained by the near detector. 


5.2. Constraints by the T2K near detector measurements 

The uncertainty of neutrino cross sections is not small, especially in the ~GeV energy 
region. Although the uncertainty is typically at the ~20 % level, the hrst order 
uncertainties of the beam and cross sections can be canceled by adopting the two detector 


technique, as described in Section [Alj For this purpose, MINOS, MINOS-I-, NOvA and 
T2K have sophisticated near detectors which collect large amounts of neutrino data to 
measure the neutrino beam flux and neutrino cross sections with high precision. 
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Figure 19. T2K 68 % (dashed) and 90 % (solid) CL regions contours of oscillation 
parameters sin^ 023 versus Am |2 for normal (top) and inverted (bottom) mass 
hierarchy |34j . The best fit point is shown as the black mark in the normal mass 
hierarchy. The Super-Kamiokande mi and MINOS [Slj results are also shown for 
comparison. 
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Figure 20. Fractional uncertainty of the T2K neutrino beam flux: (left) muon 
neutrino component and (right) electron neutrino component [^. 


CCOn-Data CClil-Data CC Other-Data 



(MeV/c) (MeV/c) (MeV/c) 


Figure 21. The muon momentum versus the scattering angle for CCOtt (left), CCItt 
( middle) and CC Other (right) samples in T2K data 


In the T2K experiment, the near detector called ND280 [13] is located at 280 m 
from the beam production target, in the same direction as the far detector, at 2.5 degrees 
off-axis. ND280 consists of two Fine Grained Detectors (FGD), three Time Projection 
Ghambers (TPG), an Electromagnetic Galorimeter system (EGAL), Side Muon Range 
Detectors (SMRD) and a 7r° detector (POD). Except for the SMRD, the detectors are 
located inside a dipole magnet of 0.2 T magnetic held. For neutrino energies around 
1 GeV, the dominant neutrino interaction is charged-current (GG) quasi-elastic (QE) 
scattering, and the second dominant one is GG 1 tt production. In the higher energy 
region, deep-inelastic scattering (DIS) becomes dominant. In ND280, the following three 
event categories have been measured: GGOvr, GGItt and GG Other where the GGOvr 
sample is for GG QE events, GGlvr for GG 1 vr production and GG Other for DIS. The 
distributions of muon momentum and scattering angle relative to the neutrino beam 
are shown in Figure M for data. The neutrino interaction models and the neutrino 
beam hux are tuned to match the observed distributions in Figure]^ After tuning, the 
uncertainties of the neutrino event rates are summarized in Table [H The uncertainties 
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Figure 22. The uncertainty of the neutrino events in the T2K far detector as a 
function of neutrino energy |34j . (Left) candidate events and (right) candidate 
events. 


of TT hadronic interactions in the far detector and the detector systematic error in the 
far detector are also shown. 


Table 2. Fractional uncertainties (%) of the number of neutrino events in the T2K far 
detector [34]. The uncertainties of cross sections are categorized into two parts: One 
is constrained by the ND280 measurement, and the other is independent of ND280. 


Sources 

Pe candidates 

z/^ candidates 

flux -|- cross sections (ND280 constrained) 

3.2 

2.7 

cross sections (ND280 independent) 

4.7 

5.0 

TT interactions in the far detector 

2.5 

3.0 

Far detector systematic 

2.7 

4.0 

Total 

6.8 

7.7 


In T2K, the number of observed v, 
uncertainty of 6.8 (7.7) 9 


(z/^) events is 28 (120) with a systematic 


In Figure 22, the uncertainties on the expected energy 
distributions of both Vf. and events are shown before and after constraint by the 
ND280 measurement. Today, the T2K sensitivity is not very limited by systematic 
errors which will also rapidly improve for the future CP violation measurement. In the 
near future, T2K expects that the total systematic uncertainty can be reduced down to 
5 % or less. 


5.3. MINOS Systematics 

MINOS measurements have signihcant statistical error. For example, in the combined 
beam and atmospheric analysis a total of 3117 beam-generated contained-vertex charged 

































Neutrino Oscillations with the MINOS, MINOS+, T2K, and NOvA Experiments. 25 

current events, distributed across the entire neutrino energy spectrum, are used lig. 
The appearance analysis hnds 172 events. Notwithstanding these small samples, MI¬ 
NOS has performed complete and detailed analysis of the systematics of the parameter 
measurements. The sensitivity of the results to systematics also benehts from the great 
similarity between the MINOS near and far detectors. Here we summarize the most 
important beam-related systematics and their effects. 

The MINOS 3-neutrino combined disappearance and z/g appearance paper isa 
describes the use of 32 systematic effects as nuisance parameters in the hnal hts. Of these 
13 concern the atmospheric neutrinos; these are not discussed further in this section. 
There are 4 dominant systematic effects in z/^ disappearance in MINOS. These are, with 
representative values (from [12]): 

(i) Hadronic shower energy (7% in the oscillation maximum region), 

(ii) /i energy (2 -3%, depending on technique), 

(iii) Relative normailzation (1.6%), and 

(iv) Residual neutral current (NO) contamination (20%) 

The systematic knowledge of the muon energy includes measurements by range 
(in the MINOS steel, 2%), and by momentum extracted from curvature in the MINOS 
magnetic held (3%). The relative normalization error of 1.6% is derived from knowledge 
of hducial masses and relative reconstruction efficiencies 021 . 

Additional systematics are taken into account for z/g appearance, affecting both 
signal and background predictions which are compared with the observed data. Many 
of the large number of systematic checks have small or negligible effects on the 
measured parameters; nevertheless they are incorporated in the hnal htting procedures 
as discussed above. 

Illustrative values of the systematic ehects are quoted in |35|, where the ehect of 
the relevant uncertainties on the far detector background prediction for z/g appearance 
in the beam are given as: 

(i) Energy scale: This includes both relative energy scale diherences between the 
near and far detectors and the absolute energy scale. The former ahects the z/g 
background prediction by directly impacting the extrapolation from data, and is 
the most important single systematic (2%). The latter enters via its ehects on 
the event selection process, and is less important. The combined ehect on the 
background is 2.7% 

(ii) Normalization: This term refers to ehects relating to the relative hducial masses 
and exposures of the two detectors. It is quoted as 1.9% 

(iii) Vr cross section: A poorly known pseudo-scalar form factor 071 causes a background 
uncertainty of 1.7%. 

(iv) All others: Small ehects due to, for example, neutrino huxes and cross-sections 
(which largely cancel due to the functional identity of the near and far detectors) 
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and hadronic shower modeling in nentrino interactions. |5H]. The snm of these 
small effects is < 1%, showing once more the effectiveness of extrapolating from 
near detector data. 

The hnal nncertainty on the backgronnd is 4%, to be compared with its 
statistical nncertainty of 8.8%. The nnmbers cited here apply to the beam mode 
only, with similar, bnt slightly higher valnes for the beam mode. In addition, there 
is a systematic error of approximately 5% on the appearance signal selection efficiency, 
stndied with CC events in which the mnon data has been replaced with a simnlated 
electron. 

6. Additional Measurements 

6.1. Additional Measurements in MINOS 

In addition to the primary mission of MINOS and MINOS+, which is the nnderstanding 
of the three-nentrino oscillation sector, the MINOS detectors and NnMI beam line are 
capable of a wide variety of additional measnrements which enrich onr nnderstanding 
of the physics of nentrinos, and other areas. Measnrements pnblished by the MINOS 
collaboration inclnde: 

(i) Searches for additional sterile nentrinos nsing charged and nentral cnrrent 
interactions m iBi, 

(ii) Measnrement of nentrino cross-sections [50] , 

(iii) Tests of fnndamental symmetries and searches for non-standard interactions. isa. 

(iv) Stndies of cosmic rays at both near and far detectors [22] • 

In this section we discnss briefly the first item, searches for sterile nentrinos, and 
present an example of a fnndamental symmetry test. 

Anomalies seen in short baseline experiments [53| and others have generated great 
interest in the possibility of a fonrth nentrino which wonld not have Standard Model 
interactions. Nentral cnrrents in the MINOS detectors are visible as hadronic showers 
withont an accompanying lepton. All active flavors of nentrinos prodnce nentral cnrrents 
eqnivalently, so that the three-nentrino oscillation phenomenon shonld not canse any 
depletion with respect to expectations in the observed far detector spectrnm. This is 
indeed seen to be the case, as docnmented in |18] which measnres a limit on the fraction 
of nentrinos which can have oscillated to sterile nentrinos, fs, of fs < 22%. More 
recently, preliminary analysis of farther data has generated limits on the sterile mixing 
angle 024 which extend the range of previons experiments (49j . 

A fnndamental test of CPT symmetry is the eqnivalence of oscillation parameters 
obtained from and The magnetized MINOS detector can perform an event-by- 
event comparison of these parameters enabling an accnrate test of this prediction. The 
resnlting allowed regions (from |16]) are displayed in fignre The difference in |Am^| 
obtained, in a two-flavor model, is |Arh^| — |Am^| = (0.12lQ2g) x 
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Figure 23. MINOS comparison of and oscillations from [15]. The 90% 
confidence level allowed regions and best fit values are shown for and oscillations, 
and for a fit in which the parameters are assumed to be identical. 


6.2. Additional Measurements in T2K 


In addition to nentrino oscillation stndies, T2K conducts various measurements on 
neutrino-nucleus cross sections. As described in Section |5.2[ the understanding of 


neutrino cross sections is important to reduce systematic uncertainties of neutrino 
oscillation measurements, which could improve the sensitivity to neutrino oscillations. 

The cross sections measured in T2K are summarized in Table As the hrst step, 
T2K measures the muon neutrino charged current (CC) inclusive cross sections with the 
T2K off-axis near detector (ND280) [T3] and the on-axis near detector (INGRID) [T7] . 
In future, these analyses will be more sophisticated to measure exclusive channels, such 
as CC-QE, CC Ivr production, and CC-coherent tt, as energy dependent differential 
cross sections. With INGRID, there are two types of neutrino detector with different 
target materials. One has an iron target, and the other has a plastic (CH) target. By 
using two targets, the CC inclusive cross sections on iron and plastic, and the ratio 
of cross sections are measured. The CC-QE cross sections are also measured with 
the CH target [TH]. With ND280, the electron neutrino CC inclusive cross sections 
can also be measured using the powerful particle identihcation performance 

of TPC and ECAL. In the analysis [20], the electron neutrino contamination in the 
beam is measured relative to the prediction in the simulation. T2K also has divided the 
measurement into two contributions: One is the electron neutrino from kaon decay and 
the other is from muon decay. The neutral current (NC) gamma production cross section 
in neutrino-oxygen interaction has also been measured PD] by using the far detector. 
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Super-Kamiokande. All the results are consistent with the predictions in the neutrino 
interaction generator libraries NEUT |2T] and GENIE [22]. Finally, with ND280, the 
electron neutrino disappearance sample was searched to investigate neutrino oscillations 
to sterile neutrinos [23] . 


Table 3. Neutrino Cross Sections measurements in T2K for charged current (CC) 
and neutral current (NC) inclusive processes. The measurements of cross sections are 
given per nucleon. The ratio of cross sections is also shown in some measurements. 


Mode 

Results (—or the ratio) 

<E^> (GeV) 

Reference 

CC inclusive 

(6.91 ± 0.13(stat) ± 0.84(syst)) x 

0.85 

m 

Ve CC inclusive 

(1.11 ± 0.09(stat) ± 0.18(syst)) x IQ-^s 

~1.3 

m 

Vfj, CC inclusive on Fe 

(1.444 ± 0.002(stat)t!(;(^?(s?/st)) x IQ-^s 

1.51 

m 

Vfi CC inclusive on CH 

U.379± 0.009Gtat)t(|;i^?(s2/st)) x IQ-^s 

1.51 

m 

Vfi CC ratio of Fe/CH 

1.047 ± 0.007(stat) ± 0.035(syst) 

1.51 

m 

CC-QE on CH 

(10.64 ±0.37(stat)t2-03(gyg^)) x 

0.93 

mi 

Vf, CC-QE on CH 

(11.95 ±0.19(stat)+i-^7Gyst)) X 10-39 

1.94 

m 

NC 7 

(1.55;°J3) X 10-33 

0.63 

m 

Ue CC ratio 

*^^'eCc/*^prediction — 1.01 dz 0.10 

~1.3 

m 

Ve (K decay) CC ratio 

*^i.'eCc/*^prediction — 0.68 dz 0.30 

— 

m 

Ve in decay) CC ratio 

*^i.'eCc/*^prediction — 1.10 dz 0.14 

— 

m 


7. The Frontier and futnre sensitivity 


7.1. NOvA 


The NOvA experiment is the principal appearance-mode long-baseline experiment at 
Fermilab. It makes use of the off-axis NuMI beam as discussed in previous sections. The 
far and near detectors minimize passive mass, moving toward the ideal of a totally active 
hducial volume which will be realized with future detectors, such as e.g. DUNE liquid 
argon TPC’s. A totally active detector typically has better performance because all 
charged particles in neutrino interactions are reconstructed properly with good efficiency. 
In the case of NOvA the active medium is liquid scintillator mixed into oil. The basic 
segmentation of the detector is planes of 3.6 cm x 5.6 cm PVC tubes, separated by 
thin walls. Light is collected in each tube by a looped wavelength shifting hber that is 
directed onto a single pixel of an avalanche photodiode (APD). The APD’s are cooled 
by a hybrid thermoelectric/water system to a temperature of —15° C. 

There are a total of 344,064 tubes arranged in an interchanging pattern of horizontal 
and vertical planes, giving an active structure 15.6 m x 15.6 m, and 60 m long. The 
active scintillator mass is 8.7 ktonnes with 5.3 ktonnes devoted to the PVC support 
structure, hber readout, and other structural components, giving a total far detector 


mass of 14.0 ktonnes. Figure shows a typical NOvA far detector event. 

In addition to the far detector, a functionally equivalent near detector is located at 
Fermilab. This detector consists of horizontal planes 4.0 m x 4.0 m, together with a 
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Figure 24. Charged-current neutrino interaction observed in the NOvA Far Detector 
(Courtesy NOvA Collaboration). 


steel muon catcher to assist in measurement of the spectrum of CC events. The expected 
average occupancy of the near detector is ~ 6 events per spill at full intensity, which 
will be separated by their time of occurrence, as is already successfully done in MINOS 
and MINERvA. 

The relatively large value of ^13 provides the opportunity for a rich harvest of 
physics results for NOvA via both disappearance and appearance, with both neutrino 
and antineutrino beams. Here we focus on the appearance channel, with its rich in¬ 
formation about 013, mass hierarchy, and dcp- Disappearance measurements are also 
highly sensitive and give excellent information on the octant of 023- During early run¬ 
ning, combination with MINOS-I- will be exploited, as discussed separately below. 

As discussed in earlier sections, the measurements of appearance of NOvA are sig- 
nihcantly affected by matter effects, potential CP violation, mass hierarchy, and the 
distinctions between neutrinos and anti-neutrinos. A useful tool for understanding 
these dependencies is the biprobability plot, in which one axis displays the appearance 
probability for i/g and the other the probability for z7g. Figure shows this situation 
graphically, for the case of maximal mixing. The results are clearly separated for many 
values of 6 cp, with areas of overlap for the regions around 6 cp = l^. { 6 cp = 37 r/ 2 ) for 
normal (inverted) mass hierarchy. Moreover, the subleading effects cause a separation 
in the appearance probabilities depending on the octant of 023. The sensitivity of NOvA 
to the octant depends on the value of 623 , and somewhat on both the hierarchy and the 
value of 6 cp. For a value of sin^2023 of 0.95, there is considerable sensitivity for all 
values of those parameters, exceeding 95% CL for signihcant portions of the parameter 
space with the nominal exposure, shown in Figure 
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Figure 25. Biprobability plot showing the effect of mass hierarchy, Scp, and vs 
exposure. The left panel shows the expected appearance probabilities for Ve (pO 
when the mixing angle 023 = r jA (maximal mixing.) The right panel shows the same 
for assumed non-maximal mixing with 023 < r/A and 023 > 7r/4. (Courtesy NOvA 
Collaboration). 


NOvA octant determination 



NOvA octant determination 



Figure 26. NOvA sensitivities for resolution of the octant of non-maximal 023, 
assumed to be at sin^ 2023 = 0.95. Left (right) panel shows the sensitivity for 023 < 7^/4 
(023 > 7r/4). (Courtesy NOvA Collaboration). 


As suggested by Figure]^ a principle goal of NOvA is to gain information about 
the mass hierarchy of the neutrino eigenstates. We can see the sensitivity of a ’’stan¬ 
dard” expected exposure in the NuMI beam in Figure ^ Hierachy and CP violation 
information are coupled, leading to the right panel of the hgure, in which the fractional 
coverage of ” CP-space” 0 — 27r at which the hierarchies can be separated is shown as a 
function of the signihcance of the separation. 


Because of parameter ambiguities, the study of CP violation in neutrino oscillations 
is particularly challenging. Full resolution of the problem may require the very large 
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NOvA hierarchy resolution 



NOvA hierarchy resolution 



Figure 27. Sensitivity of NOvA to the neutrino mass hierarchy. The left panel shows 
the sensitivity as a function of bcpjr. The right panel shows the fraction of bcp space 
covered for various degrees of statistical significance a. Both plots are for maximal 
023- Approximate one third of the bcp range gives a hierarchy determination at 95% 
C.L. (Courtesy NOvA Collaboration). 


Example NOvA contours 36x10^° POT Example NOvA contours 36x10^° POT 




Figure 28. Examples of joint NOvA sensitivity contours for sin^ 6*23 and bcp, for 
a nominal run of 36 x 10^° protons on target in the 700 kW NuMI beam. In the 
left panel @23 is maximal, and bcp is chosen for maximum separation between the 
hierarchies. The inverted hierarchy and some values of bcp are disfavored. A less- 
favorable case is shown in the right panel, where bcp is chosen to illustrate the 
difficulty of distinguishing the hierarchies. In both-cases ^13 is taken as an external 
input (Courtesy NOvA Collaboration). 


detectors of Hyper-K and DUNE, currently under discussion. However, particularly in 
favorable cases, important information can be obtained. Figure illustrates the situ¬ 
ation by plotting the simultaneous signihcance of two quantities, sin^ 6^23 and bcp, for 
both hierarchies. In confused cases information is obtained about the likely correlations 
of bcp and the hierarchy, and less ambiguous cases will favor specific regions of bcp- 


7.2. MINOS+ 

The MINOS experiment, as described in this paper, has studied the region in L/E 
near the oscillation minimum in detail, using primarily low-energy settings of the NuMI 
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Figure 29. Ratio of oscillated to unoscillated predictions for charged-current events 
in the MINOS-f experiment, as a function of exposure in terms of protons on target 
(POT). The statistical precision in the neutrino energy region of 5-7 GeV is much 
improved by the high rates available from the NuMI complex (Courtesy MINOS-f 
Collaboration). 


beam. As discussed in section 2.2 the requirements of off-axis kinematics for the NOvA 
experiment lead to a need for a higher energy on-axis setting of the on-axis beam. This 
beam, with its associated larger event rates in both MINOS near and far detectors, 
is exploited by the MINOS-b experiment. Fig shows the expected structure of the 
data spectrum which will be obtained by MINOS-I-. The experiment, will collect on the 
order of 3000 CC and 1200 NC events for each exposure of 6.0 x 10^° protons on the 
NuMI target (roughly annually). These large event rates allow a varied physics program. 


Representative physics goals of the MINOS-I- experiment include: 

(i) Precise verihcation of the expected spectral shape of the oscillation phenomenon, 
using the disappearance technique, especially in its transition region in energy 
between 5 and 7 GeV. 

(ii) Utilization of the precise spectrum together with MINOS and NOvA data to 
continue improving understanding of the oscillation parameters. 

(iii) Further study of the possibilities of sterile neutrinos, using both CC and NC 
disappearance. 

(iv) Improved precision on searches for exotic phenomena. 

Figureshows an example of (ii), considering all the NuMI program data expected 
to be obtained in the 2015 time frame. Comparison with Figure [TB] shows improved 
determination of the parameters, especially in the case of the normal mass hierarchy. 
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Figure 30. Expected near-term precision of determination of the oscillation 
parameters sin^ 023 and Am |2 for the MINOS-I- experiment. The left panel shows the 
results when MINOS and MINOS-f data are combined, and the right panel includes 
the expected data from the NOvA experiment. The contours are generated using as 
inputs the best-fit MINOS values from [54] (Courtesy MINOS-I- Collaboration). 


7.3. Future prospects of T2K 

The approved beam for the T2K experiment is 7.8 x 10^^ protons on target (POT). 
The results of T2K reported in this paper are based on 6.6 x 10^° which is only 8 % of 
the original goal In the near future, the J-PARC accelerator plans to increase the 
repetition rate of the acceleration cycle by updating the power supply system. With 
the upgrade, the beam power of J-PARC will reach 750 kW, and T2K will accumulate 
the design beam within several years. In this section, we show the physics sensitivity of 
T2K with 7.8 x 10^^ POT. In T2K, there are two beam operation modes: one is neutrino 
beam and the other is anti-neutrino beam. Since the fraction of the anti-neutrino beam 
to the neutrino is not fixed yet, we will show both possibilities. 

The current goals of T2K with 7.8 x 10^^ POT are 

• Initial measurement of CP violation in neutrinos up to a 2.5a level of significance. 

^ T2K collected 1.0 x lO^i POT on March 26th, 2015. 
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• Precision measurement of oscillation parameters in the disappearance with 
precision of 6 Am '^2 — 10 ”^ eV^ and (5sin^26'23 — 0.01 ; also determination of 623 
octant at 90 % C.L. if | 6*23 — 45°| > 4°. 

• Contribution to the determination of the mass hierarchy. 

7.3.1. Neutrino events with 7.8 x 10^^ protons on target (POT) 

Based on the analysis method in [HIIS], the expected number of and eventsj^are 
shown in Table and using the following neutrino oscillation parameters: sin^ 20 i 3 = 
0.1, sin^ 6*23 = 0.5, sin^26'i2 = 0.8704, Am ^2 = ’^•6 x 10“® eV^, |Am 32 | = 2.4 x 10“^ eV^, 
5cp = 0, and Am |2 > 0. 


Table 4. The expected number of Vf, events in T2K with 7.8 x 10^^ POT for the 
neutrino beam mode and 7.8 x 10^^ POT for the ant-neutrino [35]. 


Beam Mode 

Total 

Signal 

Signal 

Tf, -t Fe 

Beam CC 

Ue + Fe 

Beam CC 

NC 

NC 

neutrino beam 

291.5 

211.9 

2.4 

41.3 

1.4 

34.5 

anti-neutrino beam 

94.9 

11.2 

48.8 

17.2 

0.4 

17.3 


Table 5. The expected number of events in T2K with 7.8 x 10^^ POT for each 
beam operation mode |35j . 


Beam Mode 

Total 

CCQE 

CC non-QE 

CC Z/g -b Ve 

I'lxiPfi) -t i^e(Fe) 

NC 

neutrino beam 

1,493 

782 (48) 

544 (40) 

4 

75 

anti-neutrino beam 

715 

130 (263) 

151 (138) 

0.5 

33 


7.3.2. CP sensitivity 


Since the electron neutrino appearance is sensitive to CP violation, the variation 
of the number of electron neutrino events with 5 cp parameters is shown in Figure 
In maximum, we expect a 27 % change, compared to no CP violation with ^cp = 0- 
Hereafter, we assume the beam exposure to be 50 % for the neutrino beam and 50 % 
for the anti-neutrino. We also assume sin^ 29i3 = 0.10 ± 0.005 as the ultimate 613 value 
from reactor experiments. In the case of the maximum CP violation (^cp = ~90°), the 
T2K sensitivity for ^cp = -90° is shown in Figure [32] with 90 % C.L. 


This study was conducted before T2K developed the special 7r° rejection algorithm. So, the number 


of NC background is higher, compared to the results in Section 3.2 
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Reconstructed Energy (GeV) Reconstructed Energy (GeV) 


Figure 31. The reconstruction energy of expected T2K electron neutrino appearance 
events with various Jcp parameters [35) . (Left) Neutrino beam operation with 
7.8 X 10^^ POT and (Right) Anti-neutrino beam operation with 7.8 x 10^^ POT. 



sin^20^3 


Figure 32. The 90 % contour region of the (5cp versus sin^ 20i3 plane in the expected 
T2K sensitivity with 7.8 x 10^^ POT [35]. POT are assumed to be equally distributed 
to the neutrino beam mode (3.9 x 10^^ POT) and the anti-neutrino (3.9 x 10^^ POT). 
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Figure 33. T2K expected difference between the true point (Jcf, sin^ ^ 23 ) and the 
point 5cp = 0 [35] . The map of x^ difference shown in color is calculated assuming no 
systematic errors. The solid contours show the 90 % C.L. sensitivity with statistical 
error only, while the dashed contours include the 2012 T2K systematic error. (Left) 
the contour in the case of Normal mass hierarchy, and (Right) in the case of Inverted 
mass hierarchy. POT are assumed to be equally distributed to the neutrino beam 
mode (3.9 x 10^^ POT) and the anti-neutrino (3.9 x 10^^ POT). 


In reality, the sensitivity to CP violation also depends on 6 ^ 23 . In Figure [3^ as the 
T2K sensititivy, we show the difference between the true point with (^cp, sin^ ^ 23 ) 
and the hypothesis test point with ^cp = 0 . 

7.3.3. Precision measurements of neutrino oscillation parameters 


Since most of the T2K measurements are limited by statistics, more data improve 
the precision of the measurements. Among the oscillation parameters, sin^ ^23 and 
|Am| 2 | are interesting because of their relatively larger uncertainty compared to other 


parameters. Figure 34 shows the expected precision of sin^ 6^23 and |Am 32 | as a function 
of POT for the normal mass hierarchy case. Hereafter, the total exposure in T2K is 
assumed to be 7.8 x 10^^ POT which are equally distributed to the neutrino beam mode 
(3.9 X 10^^ POT) and the anti-neutrino (3.9 x 10^^ POT). The statistical uncertainty of 
sin^ 6*23 and |Am 32 | is 0.045 and 0.04 x 10“^ eV^, respectively, at the T2K full statistics. 
The precision of sin^ ^13 is influenced by the precision of other oscillation parameters 
including sin^ 6^23 and 5cp. 

An interesting question about 6^23 is which 6^23 is exactly 45° or not. In the case of 
6^23 7 ^ 45°, which octant does the value of 6^23 ht in, ( 6^23 > 45° or 6*23 < 45°)? Figure 35 
shows the region where T2K can reject the maximum mixing ^23 = 45° and the regions 
where T2K can reject one of the octants of 6 ^ 23 . The octant of 6^23 is determined at 90 % 
C.L. if 1^23 - 45°I > 4°. 
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POT POT 


Figure 34. T2K expected sensitivity to sin^ 623 and |Am§ 2 | as a function of POT for 
the normal mass hierarchy |35j . The solid lines are without systematic error and the 
dashed one with the systematic error from the T2K 2012 analysis. POT are assumed 
to be 3.9 X 10^^ POT for the neutrino beam mode and 3.9 x 10^^ for the anti-neutrino. 
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Figure 35. (Left)The shaded region is where T2K has more than a 90 % C.L. 
sensitivity to reject maximal mixing [35]. (Right) The shaded region is where T2K 
has more than a 90 % C.L. sensitivity to reject one of the octants of 823 |35|. The mass 
hierarchy is considered unknown in the normal mass hierarchy case. The shaded region 
is calculated assuming no systematic errors (statistical error only), and the dashed 
contours show the sensitivity including the systematic errors. POT are assumed to be 
3.9 X 10^^ POT for the neutrino beam mode and 3.9 x 10^^ for the anti-neutrino. 
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Figure 36. The regions where the wrong mass hierarchy is expected to be rejected at 
90 % C.L. by the NOvA measurement (blue line) and the NOvA + T2K measurements 
(black shaded regions) [3S]. (Left) the true hierarchy is normal and (right) it is inverted. 
The T2K POT are assumed to be 3.9 x 10^^ POT for the neutrino beam mode and 
3.9 X 10^^ for the anti-neutrino. The NOvA POT (3.6 x 10^^) are also assumed to be 
distributed to the neutrino beam mode and the anti-neutrino equally. 


7 . 3 . 4 . Mass Hierarchy 

Because of the relatively short baseline (~300 km) of T2K, the experiment is less 
sensitive to the mass hierarchy (more sensitive to CP). However, the measurement 
of T2K (sensitive to CP) can contribute to improving the mass hierarchy sensitivity 
of NOvA by helping untangle the two effects of CP and mass hierarchy in NOvA. 
Figure shows the 90 % sensitivity region for mass hierarchy with the T2K and NOvA 
measurements. The sensitivity is really expanded by adding the T2K measurements, 
especially for ^cp ~ 90°(—90°) case in the normal (inverted) hierarchy case. 

8. Conclusions; Building the future 

We have presented the results and prospects of neutrino oscillation measurements 
by the present generation of experiments: MINOS/MINOS+, T2K and NOvA. The 
phenomenology of neutrinos is being rapidly revealed by these experiments with a 
dramatical improvement of the precision of neutrino oscillation parameters in 10 years. 
In the standard neutrino oscillation scenario, all three mixing angles have been measured 
and found to be large enough to explore CP violation. Surprisingly, the new data from 
the accelerator experiments is beginning to become sensitive to CP violation, when 
coupled with the precise knowledge of the mixing angles. 

Upcoming results from on-going experiments, especially T2K and NOvA, will have 
large impact for the following reasons. 

• The measurements of T2K and NOvA individually are the most sensitive to CP 
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violation. By combining both, the sensitivity will be further improved. 

• In order to explore CP violation, the precision of mixing angles is essential. In 
particular, the value of 6*23 plays a key role in specifying the complicated parameter 
space of delta-CP and the mass hierarchy. These experiments are the most sensitive 
to 922 .- 

• The NOvA experiment has some sensitivity to the mass hierarchy. By combining 
with T2K, the sensitive region will be expanded. 

• The future neutrino experiment, Hyper-Kamiokande and DUNE, will have greatly 
expanded sensitivity to CP violation. The experiences of T2K and NOvA, together 
with the improvement of systematic uncertainties are key inputs to the future 
experiments. 

Going beyond the standard neutrino oscillation scenario, unexpected phenomena 
may appear in the most sensitive experiments. Thus, there is discovery potential for 
MINOS/MINOS+, T2K and NOvA at any time. 

Success of on-going experiments is building the future. 
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